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Topological matter in 3D is characterized by the presence of a topological BF term in its long- 
distance effective action. We show that, in 3D, there is another marginal term that must be added 
to the action in order to fully determine the physical content of the model. The quantum phase 
structure is governed by three parameters that drive the condensation of topological defects: the BF 
coupling, the electric permittivity and the magnetic permeability of the material. For intermediate 
levels of electric permittivity and magnetic permeability the material is a topological insulator. We 
predict, however, new states of matter when these parameters cross critical values: a topological 
superconductor when electric permittivity is increased and magnetic permeability is lowered and 
a charge confinement phase in the opposite case of low electric permittivity and high magnetic 
permeability. Synthetic topological matter may be fabricated as 3D arrays of Josephson junctions. 
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I. INTRODUCTION 

The last decade has seen a surprising development. 
Topological field theories [1], originally invented as 
background-independent theories of quantum gravity, 
have found a beautiful realization in strongly correlated 
condensed matter states. 

The general idea, pioneered by Wen [2|, is to write 
the conserved matter current describing the fluctuations 
about such topologically ordered states as the curl of a 
gauge field and to formulate the low-energy effective field 
theory for the fluctuations as a gauge theory. In two 
spatial dimensions, the dominant gauge action at long 
distances is the famed Chern-Simons term 

[3, 4], which 

describes incompressible quantum Hall fluids and their 
edge excitations. 

Based on exactly the same idea we showed 5, 6] how 
effective field theories with a topological BF term [lj pro- 
vide a unified description of topological matter in any 
space-time dimensions. The relevance of the BF term for 
topological insulators has been recently derived in 0, [1| . 

The BF term [l[ is a P- and T-conserving, higher- 
dimensional generalization of the Chern-Simons terms. 
In two spatial dimensions it reduces to a mixed Chern- 
Simons term, which, for non-compact gauge groups, can 
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be diagonalized into the sum of two standard Chern- 
Simons terms of opposite chiralities: such theories have 
often been referred to as doubled Chern-Simons mod- 
els in the literature 0, 0, 0). In two spatial dimensions, 
the topological Chern-Simons term dominates the long- 
distance effective theory. In three spatial dimensions, in- 
stead, there is another marginal term that co-determines 
the infrared behavior and must be added to the action: 
this is the standard Maxwell term. 

What distinguishes the different possible phases of this 
model is the role of the topological defects ensuing from 
the compactness of the U(l) gauge group. The compact- 
ness of this group and its consequences are often over- 
looked in the literature on this subject. It requires an 
unavoidable UV cutoff in the formulation of the model 
10] and it is this UV cutoff that prevents the diagonal- 
ization of compact mixed Chern-Simons terms Q. 

In this paper we show that the T — quantum phase 
structure of topological matter is determined by three 
parameters: the electric permittivity and the magnetic 
permeability of the material and the topological coupling 
constant. For intermediate levels of electric permittivity 
and magnetic permeability the material is a topological 
insulator [7[. We predict, however, new states of mat- 
ter when these parameters cross critical values: a topo- 
logical superconductor when electric permittivity is in- 
creased and magnetic permeability is lowered and a dual 
charge confinement phase in the opposite case of low elec- 
tric permittivity and high magnetic permeability. 

Our results raise the possibility of turning a topological 
insulator into a superconductor, if a mechanism can be 
devised to lower the magnetic permeability and increase 
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the electric permittivity, while preserving the topological 
entanglement of the ground state. Equally interesting 
would be the realization of a charge confinement phase 
for high magnetic permeability and low electric permit- 
tivity. We also point out that these results can be studied 
in synthetic BF matter, constituted by fabricated three- 
dimensional Josephson junction arrays. 

Let us write the conserved charge fluctuations above 
the topological ground state in terms of a Kalb-Ramond 
[HI two-form gauge field bp V as j M = (k/2Tr)e liUa pd v b o , . 
We consider also an effective statistical gauge field 
describing the braiding of charged particles and vortex 
lines. The corresponding field strength is the pseudo- 
tensor </> M „ = {k/2ir)€p V apdaap describing the current of 
vortex excitations. Then the BF term is the topological 
action J d 4 x aptp Va0 d v b a ^. Since and bp V have 
canonical dimensions [mass] and [mass] 2 (we use units 
c = h = 1), respectively, k is the dimensionless BF cou- 
pling. In applications to topological matter the BF term 
represents the effective action for the matter degrees of 
freedom, obtained by integrating out all microscopic in- 
teractions, including electromagnetism. 

Contrary to the case of two spatial dimensions, how- 
ever, in three spatial dimensions, the topological BF term 
is not the only action term dominating the infrared be- 
havior. Another marginal term can be added to the ac- 
tion: this is a Maxwell-like term constructed from the 
effective gauge field a M . We shall thus consider con- 
densed matter systems with the following long-distance 
Euclidean effective action: 



5' 



ik 
2^ 



d i xa fl e lxua pd v b a 0+ / d 4 x 



1 



1 



2e 2 A 



2e 2 T] 



(1) 

where e and b are the effective electric and magnetic 
fields constructed from the gauge field a M , e is the elec- 
tron charge and r\ and A are two dimensionless couplings. 
The effects of possible topological term cx J d 4 xe ■ b will 
be discussed elsewhere. 

In what follows we will also add to the action an 
infrared- irrelevant kinetic term for bp V : 
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with = \ep Va pf va p, fp Up = dpb vp + d v b pfi + d p bp V the 
Kalb-Ramond field strength, and g a parameter with the 
dimensions of a mass. This serves only as an ultraviolet 
regulator for Gaussian integrals. It describes the massive 
modes that live on short length scales <C — = ■; — ZL ^=. 

First of all let us analyze the meaning of the pa- 
rameters A and r/. To this end we introduce the cou- 
pling with an external electromagnetic field: ijpA p = 
T^A^e^vapdvbap, in (JTJ Gaussian integration over a p and 
b a p induces the electromagnetic effective action. At long 
distances (^> -M this is given by: 
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where E is the physical electric field and B the magnetic 
one. This shows that the model (JTJ describes a topo- 
logical insulator and identifies the parameters of the low- 
energy effective action as the inverse electric permittivity 
and magnetic permeability, respectively, 
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One of the main characteristics of topological in- 
sulators is the appearance of the so-called axion- 
electrodynamics term [12, [r| : 
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in their action when time-reversal symmetry is broken, 
the parameter 9 being quantized in odd multiples of n: 
9 = (2n + 1)tt ,n £ Z [12]. Fp V is the field strength of 
Ap and F^ v its dual. 

We will now show how this term arises naturally in our 
model due to a (T-symmetry breaking) direct coupling of 
the statistical gauge field to electromagnetic fields. This 
is represented by the additional term -gt^4>p V Fp U in the 
action ([1]). Gaussian integration over the matter fields 
dp and b^v induces, at long distances (3> — ), the axion 
electromagnetic term with parameter = 4^. Requiring 
a T- invariant partition function exp(— S) is equivalent to 
the Dirac quantization condition M = integer. Note that 
the BF coupling k plays the same role as the charge of 
the condensate in spontaneous symmetry breaking. As 
we show below k is indeed an integer. 

The action ([T]) has two U(l) gauge symmetries under 
the transformations: 
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The important point is that these gauge symmetries are 
compact and the compactness of the U(l) gauge group, 
as usual, leads to the presence of topological defects. As 
is beautifully explained in [lfj, the compactness of U(l) 
gauge groups inevitably introduces a length scale, the 
compactification radius of the gauge fields. There are 
essentially only two ways to introduce this UV cutoff in 
a gauge invariant manner: either one embeds the U(X) 
group in a larger, non-Abelian compact group and one 
breaks the symmetry down to U(l) on a given scale or 
one formulates the entire theory on a lattice. We choose 
the latter. 

The lattice we consider is a hypercubic lattice in four 
Euclidean dimensions, with lattice spacings I. On this 
lattice we define the following forward and backward 
derivatives and shift operators: 



= f(x+lp,)-f(x) 
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Spf(x) = f{x + lfi) , 
Spf(x) ^ f{x - IpL) . (7) 



Summation by parts interchanges the two derivatives, 
with a minus sign, and the two shift operators. We also 
introduce the three-index lattice operators [5j: 
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These operators are gauge-invariant in the sense that: 
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Moreover they satisfy the equations: 
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surfaces. Its partition function can be easily obtained by 
a Gaussian integration over the gauge fields and bp V 
in PJ) : 

^Top = X) fQ^} exp( — ^Top) , 

{M M1/ } 

+<? 2 fc 2 A Moi^^Moi + g 2 fc 2 ?7 Af J3 ^M, 3 

+^ Q wft^ Map , (13) 
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The expressions Op^xu and Dp V are lattice versions of 
the Kalb-Ramond and Maxwell kernels, respectively. 

On the lattice, the gauge fields a M and bp V become link 
and plaquette variables, respectively. The compactness 
of the gauge groups implies that they are angular vari- 
ables, invariant under the shifts: 
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a.p -> a M + — n M , bp V bp V + -pUp V ; rip,, Up V G Z , 

(11) 

where =f- is the radius of the U(l) group. In order to im- 
plement this invariance one introduces in the Euclidean 
partition function integer link and plaquette variables Qp 
and Mp V as follows: 

z = 12 {«*} / Va v I Dbfuj, exp(-S) , 

+ 2g~ T X~f i f i ~ i^2n a H^Hapb a p + 

+ik(lQpap) + ik{lMp V bp V ) , (12) 

where fp V = Kp Va a a and fp = \Kp aj3 b a p . 

The integer-valued variables Qp and Mp U appearing 
in (p~2]) represent topological excitations whose role is to 
make the BF term periodic. They satisfy dpQp — and 
dpMp V — dvMp V = 0. In Euclidean space they represent 
closed electric lines and compact magnetic surfaces. The 
corresponding physical excitations in Minkowski space 
are the world-lines of quasiparticle charges Qg and the 
world-surfaces of magnetic vortex lines M oi . These topo- 
logical excitations are the singularities in the two field 
strengths due to the compactness of the two gauge sym- 
metries. It can be shown [l4| that they represent char- 
geons and spinons in 3D. As we now show, the physics of 
a compact, U(l) BF term is much richer than that of its 
non-compact cousin with gauge group K. 

The phase structure of our model is determined by the 
statistical mechanics of a coupled gas of lattice loops and 
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wth the mass m defined by m 2 

(X/r))dodo + didi the lattice Laplacian. 
l/y/ejl is the speed of light in the topological medium. 
This partition function can be interpreted as the Eu- 
clidean partition function for a lattice model of interact- 
ing particles (whose world-lines are parametrized by the 
closed loops Qp) and closed vortex strings (whose world- 
sheets are parametrized by the compact surfaces Mp U ). 
By choosing the analogue of the Landau Ginzburg pa- 
rameter a = ml > 0(1) the Yukawa interactions in (TIB"]) 
reduce essentially to contact terms expressing the self- 
interactions of topological defects. If the imaginary terms 
in the action are absent, the model (fT3"j) can be treated 
by standard statistical mechanics techniques. This is re- 
alized when k its an integer, k — n, n G Z. 

In order to derive the quantum phase structure of the 
model at T = we consider the free energy of static 
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0) charges and vortex lines. This is given 
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where Sq and Sm denote the entropies of closed lines 
and compact surfaces on the lattice, respectively.These 
are proportional to the length of world-lines (for quasi- 
particles) and assumed as proportional to the area of 
world-surfaces for vortex lines. The exact form of their 
entropy is not known although these are the well-defined 
confining strings introduced by Polyakov [l5|. We can 
thus derive only qualitative results. 

For k 2 X small, the entropy of electric strings dominates 
their self-energy and the electric part of the free energy is 
minimized by configuration in which long electric strings 
condense. The contrary is true for large k 2 X in which case 
electric strings are dilute. The same reasoning indicates 
the condensation of magnetic surfaces for large values of 
77, whereas they are dilute for small r\. 

All together we find an electric condensation phase for 
k 2 X <C a 2 /e 2 and r\ <C 7r 2 /e 2 and a magnetic condensa- 
tion phase for k 2 X ^> a 2 /e 2 and 77 7r 2 /e 2 . In between 
we have the original topological insulator phase, in which 
both types of topological excitations are dilute. 

In order to distinguish the various phases we introduce 
the typical order parameters of lattice gauge theories [T3| 
namely the Wilson loop for an electric charge q and the 
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't Hooft surface for a vortex line with flux 



Lw = exp iq 



HI 



Sh — exp i(, 



where j° x vanishes everywhere but on the links of a closed 
loop, where it takes the value 1, and where (jf* v vanishes 
everywhere but on the plaquettes of a compact surface, 
where it takes the value 1. Since the loops are closed and 
the surfaces compact they satisfy d^j^ — d^eff^ = 0. 

The expectation values of (Lw) and (Sh ) can be used 
to characterize the various phases. These expectation 
values are given by 



(L w ) 
(Sh) 



^To p (Q^ + | 



Z T o P (Q fl ,M t ,„) 



(16) 



where the notation is self-explanatory. The dominant 
contributions are given by Q 



perfect spinon conductivity, electric screening and a pho- 
ton mass fl(| in the magnetic condensation phase. We 
thus conclude that the electric condensation phase is a 
topological superconductor, whereas the magnetic con- 
densation phase is a charge confinement phase (spinon 
or dual superconductor). As explained above these topo- 
logical superconductors are characterized by a flux quan- 
tization condition 4> = 2i:/k. The charge confinement 
phase is the exact dual of the superconducting phase 
|17| . The photon becomes massive due to the Stiickelberg 
mechanism, electric fields are expelled and restricted to 
thin flux tubes between particle-antiparticles pairs and 
no charged excitations can be observed. Note that the 
original 27r periodicity in 9 of topological insulators is 
lost in the confinement phase. All together this amounts 
to the following T = quantum phase structure 



intermediate regime 
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top. insulator ; 

top. confinement . (19) 
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As is evident from these expressions, the effect of the 
condensation of topological defects is to transform the 
original short-range Yukawa interactions of chargeons 
and spinons into long-range Coulomb interactions. The 
consequences of this are best appreciated by express- 
ing the 't Hooft surface and the Wilson loop in terms 
of external gauge potentials and B^, respectively, 
</>°* oc K^ va A a and j® x cx K^ufiBap. Then the induced 
charge and vortex currents can be computed as: 
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These are London equations for the induced charge and 
spin currents, which express perfect charge conductivity 
and a photon mass in the electric condensation phase and 



As we have shown in Q, topological matter can be 
simulated in synthetic materials constructed as arrays 
of Josephson junctions. In what we called the self-dual 
approximation, such arrays are governed exactly by the 
model (O with k = 2 and: 
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where Ej and Ec are the Josephson coupling and the 
charging energy of the array, respectively. In this case 
the phases correspond to the superconducting, metallic 
and insulating phases of the array. The phase transi- 
tions have to be analyzed by varying Ej and Ec at fixed 
ml = ^/8EcEjl 2 > 0(1) and are due to a competition 
between the superconducting junction favoring a global 
phase lock-in and the charging energy favoring the local- 
ization of charges on the islands. 

Note also that a phase structure analogous to (fT9]) has 
been found in [l8| in the compact BF formulation of 2D 
doped Mott insulators. 
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